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Force-constant and positional disorder have been introduced into diamond lattice models in an attempt to 
mimic the vibrational properties of a realistic amorphous silicon model. Neither type of disorder is sufficient 
on its own to mimic the realistic model. By comparing the spectral densities of these models, it is shown that 
a combination of both disorders is a better representation, but still not completely satisfactory. Topological 
disorder in these models was investigated by renumbering the atoms and examining the dynamical matrix 
graphically. The dynamical matrix of the realistic model is similar to that of a positionally-disordered lattice 
model, implying that the short-range order in both systems is similar. 



1 Introduction The physics of atomic vibrations in perfect crystals is well understood |T)l2l. In the 
presence of disorder (e.g. noncrystallinity), many extra phenomena are observed: for example, localization 
1 3 1 , the boson peak (extra modes at low frequencies) 1 4 1 , and anomalies in the specific heat 1 5 1 and thermal 
conductivity |6J. The microscopic origin of these and other, similar, phenomena have been difficult to de- 
duce. In particular, it would be helpful to know the contributions due to different types of disorder. In this 
paper, we consider the effects of three types of disorder: force-constant disorder, in which the strengths of 
the springs between the atoms are disordered, positional disorder, in which the atomic positions are disor- 
dered, and topological disorder, in which the network of bonds is disordered. We compare the behaviour of 
models with force-constant and positional disorders (separately and in combination) to those of a realistic 
atomic model. A further aim is to see if it is possible to mimic the behaviour of amorphous systems with 
an appropriately-disordered crystal. 

2 Force-constant and positional disorder In this paper, we study 3D diamond-cubic crystal lattices, in 
which all atoms are connected by springs, and have the same mass. We assume that the displacements from 
equilibrium are small, and use the harmonic approximation. Hence, we can treat the system classically (fl, 
and describe the vibrations by a Hermitian dynamical matrix D. The problem can then be treated in the 
Hamiltonian formalism |7 8 , 9|, with energy e = uo"^, the squared vibrational frequency. We study the 

effect of disorder on the spectral density: A(e) = (^5{e — D)^, where (. . .) represents configurational 

averaging. We can represent the spectral-density operator in the crystalline basis |k, (3), where k is the 
wavevector of a plane wave, and (3 is the polarization 1 10|. In this case, the spectral density A\^p{e) is: 




Ak/3(e) = ( > J(d|k,/3)|^<5(e-e,) ), (1) 

where | ((i|k, /3) p are the weights of the (disordered) eigenstate \d) (with energy Cd) in the crystalline 
eigenstate |k, /?) . Large (~ 10^ atoms) diamond-cubic crystals have been simulated, with nearest-neighbour 
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Fig. 1 Longitudinal (circles) and transverse (crosses) dispersion relation for a positionally-disordered di- 
amond lattice with a = O.lri, where ri is the nearest-neighbour distance. Mixing between polarisations is 
clear. 



spring constant ki and next-nearest-neighbour spring constant K2- A ratio of ki/k2 — 8 was chosen, as 
this gave a vibrational density of states (VDOS) most like that of a crystal with the more realistic modified 
Stillinger- Weber (mSW) potential fl 1' 121. Force-constant disorder was added by randomly taking the 
spring constants from a uniform distribution of half width A, centred on the crystalline spring-constant 
value. The width A is expressed as a fraction of the crystalline spring-constant value k, and this fraction 
is the same for both spring constants. To avoid mechanical instabilities, widths were limited such that all 
force constants were positive, i.e. A < k. Positional disorder was introduced by randomly shifting all 
atoms by distances, chosen from Gaussian distributions centred on zero with half width a, independently 
in the x-, y- and z-directions. 

Figure^shows the dispersion relation for a positionally-disordered diamond lattice with a = O.lri, 
where ri is the nearest-neighbour distance. At each value of k, the spectral density has been orientationally 
averaged, and the positions of the peaks are shown in the figure. The dispersion relation shows many of 
the expected features: three acoustic branches are observed, of which one is longitudinal, and two are 
transverse. Mixing between the polarisations is seen: the longitudinal spectral density picks up the peaks 
due to the transverse acoustic (TA) branches. The optic modes are less clearly defined; the peaks due to 
these modes were broadened considerably by even small amounts of disorder Nevertheless, the presence 
of optic branches at high energy in both the transverse and longitudinal spectral densities is observed. 

The effect of a low degree of positional disorder on the spectral density is qualitatively very similar to 
the effect of small force-constant disorder 1 13); the peak shifts in position, and broadens. When compared 
to the spectral density of a realistic 4096-atom model of amorphous silicon, generated with a modified 
WWW procedure fl4' '151 and relaxed under the mSW potential, it was observed that neither type of 
disorder was adequate on its own to reproduce the features of the amorphous model's spectral density. 

A combination of these two types of disorder is better able to reproduce the spectral density, as seen 
in Figure 12] The disordered lattice models are able to mimic the main features of the amorphous model, 
with peak positions and shapes at least roughly correct, but they often have more peaks (particularly at 
high k) than the rather smoother amorphous spectral densities. This is due to the orientational averaging 
- each of these peaks comes from a different orientation of the k-vector The fact that these peaks are at 
different energies suggests that, despite the added disorder, the lattice models are still anisotropic, whereas 
the amorphous model is isotropic (since the multiple peaks are absent). This implies that a form of disorder 
needs to be considered that will render the lattice models less anisotropic. 
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Fig. 2 Orientationally-averaged spectral densities of a realistic 4096-atom amorphous silicon model com- 
pared to disordered 125 000-site diamond lattices. Each graph shows the spectral densities of the realistic 
model (solid line), a lattice model with force-constant disorder A = 0.75ft (dashed line), a lattice model 
with positional disorder u = O.lSri (dotted line), and a model with both the above force-constant and posi- 
tional disorders (dot-dashed line), k — 0.40lA^^ (fc = 0.41lA^^ for the realistic model): (a) longitudinal, 
(b) transverse; k = 1.200A~^ (fc = 1.197A~^ for the realistic model): (c) longitudinal, (d) transverse. 
The value of k is slightly different between the lattice and realistic models due to the necessity of using 
k-vectors that fit the periodic boundary conditions. 



3 Topological disorder Another type of disorder is topological disorder, in which the bond network is 

disordered such that the ring-size distribution is no longer crystalline. Topological disorder was examined 
by studying the dynamical matrix. If two atoms i and j interact, then the dynamical matrix elements 
representing that interaction will be non-zero. We make use of a simple atom renumbering scheme and 
show the dynamical matrix graphically, to compare the topological disorder present in the reaUstic model 
and in positionally-disordered models. 

Atoms were renumbered by first dividing the models into small cubes, such that the average occupation 
was one atom per small cube. The models are cubic, with side L. The origin was taken at one corner of the 
model (say x,y,z = 0), and atoms were numbered by moving up through all the small cubes above it in the 
positive ^;-direction, numbering the atoms in order of the ^-coordinate of their position. Once the z = L 
face of the model was reached, the scheme returned to the z = face, and moved along one small cube 
distance in the y-direction. As before, atoms were numbered by moving up in the ^;-direction. This process 
was repeated until the y = L face was reached, when the scheme returned to y,z = and shifted along 
one small cube in the a;-direction, and then repeated the numbering above, until all atoms were numbered. 
The atom number now corresponds (in some way) to its position. 
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Fig. 3 Graphical representation of the dynamical matrix of: (a) a diamond crystal with nearest- and next- 
nearest-neighbour interactions; (b) the same crystal positionally disordered with a — 0.6ri; (c) a realistic 
4096-atom model of amorphous silicon. 



In this renumbering scheme, the dynamical matrix of a crystalline diamond lattice with nearest- and 
next-nearest neighbour interactions is shown in Fig. [3^. The effect of introducing positional disorder 
(fj = 0.6ri), into the diamond lattice, is shown in Fig.|3j5. The realistic model of amorphous silicon, with 
the same renumbering scheme, is shown in Fig. |3l;- Remarkably, the dynamical matrix for the realistic 
model (Fig. |3j;) looks very similar to that of the positionally-disordered crystal (Fig. |3j3). We interpret 
this as implying that the short-range structural order, responsible for the vibrational behaviour, in the 
topologically-disordered amorphous model is very close to that in a positionally-disordered crystal. 

4 Conclusion Examining the spectral densities of diamond lattice models with added force-constant 
and positional disorders shows that these disorders are insufficient (either by themselves or in combination) 
to mimic completely the vibrational behaviour of a realistic model of amorphous silicon. Representing the 
neighbour list graphically shows that the short-range order, responsible for the vibrational behaviour, in 
a positionally-disordered diamond lattice is similar to that of a realistic model, and this is an avenue for 
further research. 
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